A Cosmological Tale of Two Varying Constants 
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We formulate a simple extension of general relativity which incorporates space-time variations 
in the Newtonian gravitation 'constant', G, and the fine structure 'constant', a, which generalises 
Brans-Dicke theory and our theory of varying a. We determine the behaviour of Friedmann universes 
in this theory. In the radiation and dust-dominated eras otG approaches a constant value and the 
rate of variation of a is equal to the magnitude of the rate of variation in G. The expansion dynamics 
of the universe are dominated by the variation of G but the variation of G has significant effects upon 
the time variation of a. Time variations in a are extinguished by the domination of the expansion 
by spatial curvature or quintessence fields, as in the case with no G variation. 
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There have been many studies of the cosmological consequences of allowing some of the traditional constants of 
Nature to change. These include evaluations of the effects of altering the observed value of a constant to another 
. constant value and studies of the time-evolution of 'constants' in generalisations of the general theory of relativity 
f* — ' that allow them to become space-time variables. The most studied case is that of varying the Newtonian gravitation 
I/"") , constant, G, through the Brans-Dicke (BD) scalar-tensor theory of gravity 0. Recently, following Bekenstein, [0, 
Sandvik, Barrow and Magueijo ||, have developed a theory (BSBM) which describes the space-time variation of 
the fine structure constant. In these, and other, studies of varying constants only a single constant is allowed to vary 
at one time. However, since we have no understanding of why the constants of Nature take the values that they 
do, whether they are logically independent, or even whether they all are truly constant, this restriction is somewhat 
artificial. Motivated by recent observational evidence for a time evolution of the fine structure 'constant', a, at 
*-£h ' redshifts z ~ 1 — 3.5, 5 @, we have unified the BD and BSBM theories to produce an exact theory which 
describes the simultaneous variation of a and G. This type of model also provides a framework within which to 
consider the consequences of changes in the scale of extra dimensions of space on apparent three-dimensional coupling 
f~) | constants. 

In section 2 we set up the theory and evolution equations for Friedmann universes in a theory that generalises 
general relativity to include varying a and G. In section 3 we show how to find the cosmological solutions during the 
dust-dominated eras. We find an exact solution where aG is constant during the dust era while a and G _1 both 
increase with time. We then determine analytically the coupled evolution of a and G during the radiation, curvature, 
and vacuum-energy dominated eras of cosmological expansion. From here we go on to check the solutions numerically 
and we show how in Universes like our own, with actual initial values for a and G the asymptotic behaviour is never 
reached. Instead we find constant a and G in the radiation era, slow growth of a and slow decrease in G in the dust 
epoch, constant values for both in curvature dominated universe, and constant a and decreasing G in A dominated 
epoch. Generally we find that the overall evolution of the expansion scale factor of the universe is dictated by the G 
variation and assumes the form found in the Brans-Dicke theory to a very good approximation irrespective of the a 
variation. The evolution of a is influenced by the G variation but does not differ much from that found in the BSBM 
cosmologies where only a varies. 



II. FIELD EQUATIONS 

We introduce the structure of the BSBM theory for varying a as another matter field in Brans-Dicke theory. The 
resulting theories has two scalar fields: the BD field cf> propagating variations in G, and the field if) propagating 
variations in a. The action for this theory becomes 

S = J d x^J-g I Rcj) + — C - ujbd — t — J (!) 
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where 



C = C m + £ eTO exp(-2'0) (2) 

and 

LO 

A* = --^^ ■ 

The field equations for the theory, specialised to the case of a homogeneous and isotropic Friedmann space-time 
metric containing dust and radiation perfect fluids are: 



3 ^ = ¥ + ICI exp(-2^)) + p r exp(-2^) + ^) - 3^ + " 4 ( 3 ) 
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^ + 3-0=—^ (Pm-2p^) (4) 

$ + 3Hi/j = -- exp(-2V>)Cpm (5) 
Pro + 3iip m = (6) 
pV + 40> r = 2ipp r (7) 

where p^, = ^ip 2 is the kinetic energy density for the ip fluid, with to the coupling setting the relevant energy scale 
for the "0-field. £ is defined as the ratio C em / p m averaged over all types of matter in the universe. The fine structure 
'constant' is given by (U — c = 1) 

a = a<j exp(2-0), (8) 

where ao is the present day value of the fine structure 'constant'. The present-day value of G is set equal to unity. 

We shall confine our attention to the case with £ < where the magnetic energy dominates the electric field 
energy of the matter coupling to electric charge in the universe. This places particular constraints upon the nature 
of the cold dark matter dominating the universe today. From our earlier studies, ||- Q, we know that this case 
provides a slow variation with a increasing logarithmically in time during the dust era but staying constant during 
any subsequent curvature or cosmological constant dominated era. Also, in a universe with a matter-radiation balance 
like our own, a remains constant during the radiation era except close to the initial singularity. Negative £ models are 
well behaved and correspond to the dark matter in the universe being dominated by magnetic coupling, (for example 
superconducting cosmic strings contribute £ = —1). The expansion scale factor evolution is not affected by variations 
in a to leading order. By contrast, the choice £ > creates major changes to cosmological evolution. It does not lead 
to slow increase of a with time during the dust era, as observations suggest, and the evolution of the expansion scale 
factor is affected to leading order (see for example refs. [p|JlO|] who discuss related theories for the variation of a with 
£ > and hence a < cosmological behaviour in the dust era in contrast to our discussions in ||- M and below). In 
what follows we shall investigate how the £ < evolution of the fine structure constant couples to variation of G in 
the Brans-Dicke theory. 

The constant lubd is the Brans-Dicke parameter and lj is the analogous parameter for the coupling of the ip field 
driving variations in a. We have used the facts that dust is pressureless, p — p r /S for a sea of radiation and p = p^ 
for a fluid with kinetic energy only. Equation (H) can be recast for numerical solutions 



a _ l<j> 1 
a ~ ~2d> 2 
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and eqn. (|7|) integrates to give p 7 exp(— 2if)) cx a~ 4 . Note that this is the combination that appears in the generalised 
Friedmann equation, (||). In ref. Q we showed how to deduce the solutions of these equations when G is constant. 
Here, we will extend this analysis to the new situation where both a and G vary in time. 



III. DUST ERA EVOLUTION 



From our study of the Friedmann models in BSBM theory we know that, to a very good approximation, the a 
variations do not significantly affect the evolution of the expansion scale factor a(t). The effects of varying G in 
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Brans-Dicke theories is different. No matter how slow the variation in G, a correction will occurs to the power of the 
time-variation of the expansion scale factor. In the dust era we assume the asymptotic solution for the Brans-Dicke 
(BD) flat dust model holds to high accuracy. This is an exact solution of (0) for £ = k = tjj — p 7 = and is the 
late-time attractor of the general flat BD dust solution (see refs. [|ll|Jl2| 1, p3| ) which differs only as t — > 0, where 
the solution becomes dominated by the kinetic energy of the <f> field and approaches the BD vacuum solution. Thus, 
to leading order the expansion dynamics and <fi evolution are described at late times by the exact Brans-Dicke dust 
solution with k = 0: 

a(t)cx^ 2 ^")/ 3 ;0 = o f l (10) 

p = Ma~ 3 M = (11) 

8ttM 
n(3 + 2lo B d) 

where n is related to the Brans-Dicke parameter by 

n = 2/(4 + 3uj BD ), (13) 

and M is the present density of the universe in Planck units, M ~ 10~ 123 . 

What is the asymptotic solution for a during the dust era? The relevant equation is (||) , which can now be rewritten 

as 



A(t 2 ->)=iVexp(-2^) (14) 



dt 
where 

N = -^ Pm a 3 = > N = 0, 

and —Q/uj ~ 10~ 4 is the best fit of this parameter ratio to the observations of Webb et. al. JtJPT. 
Unlike the case with constant G, there is an exact solution (for lubd positive and finite) 

^) = |ln(i) + iln7V-iln(|-y) (15) 



so we have, using this solution for ip to solve for in (||): 



2Nt n 

a(t) = a exp(2^) = a — -, r (16) 

nil — n) 

Hence, there is a simple relationship between a(t) and G(t): 

< f )=G -l = 27rtJ ( 1 ~ n ) a 2ttuj{2 + 3u> B p) a ^ f » ^ 

-C(3 + 2w B L))ao — C(3 + 2wbd)(4 + 3lo B d) a 

so aG is always a constant. Note that for large values of ujbd we have a simple relation between the values of G and 
a: 

G— » > (18) 

As expected, a increases whilst G falls as t — > oo in a flat universe. It is interesting to note that the asymptotic value 
of Ga is uniquely determined by the parameters in the model with no arbitrary constants. 

Although the asymptotic behaviour is now determined, the question of whether this can be reached on a cosmological 
timescale depends strongly on the choice of initial conditions and needs to be investigated numerically. We can quickly 
conclude that the asymptotic regime is not reached in our universe. Presently we have a ~ 1/137, and in our units the 



1 The value used for —(,/lu is the value fitted for the BSBM theory with constant G. However, since n is so close to zero it 
should not be significantly different numerically in the case with varying G 
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numerical value of M is extremely small, ~ 10 -123 . Obviously the actual value of a is then many orders of magnitude 
larger than given by the solution in eq.((l^) and we are thus nowhere near the asymptotic regime. Consequently, in 
order to find the behaviour of a and G we turn to numerical solutions of the equations. We evolve the Friedmann 
equations through time with initial conditions chosen so as to yield the present day values of G and a. For a we find a 
behaviour very similar to the BSBM theory, with a slow growth giving a relative change of the order 10" 4 throughout 
the dust epoch. G goes through a decrease of order 10 -3 during the same period. The numerical results are shown in 
Figures [l] and ||. 
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FIG. 1. Evolution of the relative shift in the values of the two 'constants' in a realistic flat cosmology with vacuum energy 
and approximately accurate initial values for the fields. We start from a radiation-dominated universe where both a and G 
stay constant. Thereafter we move into dust domination where a changes slowly, while G goes through a small decrease. As 
the universe becomes dominated by the vacuum energy, a goes to a constant, while G goes on decreasing indefinitely as in 
ordinary Brans Dicke theory. Values used for the couplings are the minimum allowed value of ujbd = 3500 and we take the 
best fit value of £/w = -10" 4 from BSBM theory. 



IV. RADIATION ERA EVOLUTION 



The evolution in the radiation era is slightly more complicated because of the contribution of the term to the 
right-hand side of the <f> evolution equation. This means that we do not have the usual late-time asymptotic behaviour 
of constant 4> to accompany the a — t x l 2 scale factor as in BD radiation universes. If we assume 



then we have 



where 



t l/2 



4(^ 3/2 ) = R-Xt 3 ^iP 2 (19) 
at 



A^-ff— , (20) 
3 + 2uj bd 

8vrM AttluN . N 

= (21) 

~ 3 + 2lj bd ~ C(3 + 2w BD ) 

The R term is negligible when the kinetic energy of the ip field dominates the matter density during the radiation 
era. Likewise, the A term can be neglected when the matter density dominates the tp kinetic energy. We also have 
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FIG. 2. Evolution of the relative shift in the values of the two 'constants' in an open universe through radiation, dust and 
curvature-dominated epochs. Initial values for the fields are set so as to give realistic values at present time. Again a and G 
are constants in the radiation dominated era, whilst a increases and G decreases through matter domination. As curvature 
starts to take over the expansion, both a and G tends to constants. Values for the couplings are the minimum allowed value 
of lobd = 3500 and the best fit value of (,/uj = — 10 -4 from the BSBM theory. 



dt 



3/2n 



iVexp[-2^] 



as in the case with constant G. This has the exact solution 

1 



i/> = - ln(8./V) + - In(i) 



as before, so ip 2 — (16t 2 ) 1 . If we substitute this in (19) 

d 



dt 



{4>t 



3/2x 



R 



,-1/2 



16 



(22) 



(23) 



(24) 



2Rt 1 l 2 -^\n{t) + Ct- 1 ' 2 
8 



(25) 



where C and </>o are constants. If the universe expands for long enough to reach the asymptotic regime then we have 
(as R > 0) 



~ ln(t) 



G- 1 



8nuiN 



so, from eqns. ( |18|) and (g), 



G— 

a 



((3 + 2ljbd) 
C(3 + 2lo bd ) - 



f l/2 



(0J B D 



8itujN 



AttluN 



(26) 
(27) 

(28) 



for large ujbd- Thus we still have the nice asymptotic behaviour of ocG in the radiation-dominated epoch. However, 
we again need to compare with numerical results to determine whether these asymptotic solutions can indeed be 
realised in the Universe. As in the case of dust, the same simple reality check can now be performed on the solution 
(p3|). As in the case of constant G we are nowhere near this particular solution in our Universe, a would need to be 
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several orders of magnitude smaller if it was to satisfy the solution, and as in the BSBM theory we expect instead 
a constant value of a in the rad epoch. This assumption is indeed confirmed by the numerical solutions shown in 
Figures (§) and (Q). 

Another possible problem for the analytic solutions above would arise if the kinetic energy of the tp field dominates 
the matter density during radiation domination. We regard this situation as unrealistic and it cannot be realised 
asymptotically. 



V. CURVATURE ERA EVOLUTION 



During a curvature-dominated phase of an open universe the expansion scale factor tends to that of the Milne 
vacuum universe, which is an exact solution of general relativity and of Brans-Dicke theory (with constant 4>) with 

a = t (29) 

Using this in the propagation equations for <fi and ip,we find that, as t — > oo, so leading order 

iVexp[-2V> ] 



ip = ip* - 



t 

A-kujN 



C(3 + 2u BD )t 



with tpif and constants, so both a and G tend to constant values as t — > oo. In a universe that passes directly from 
dust domination to curvature domination these constant values will be very close to the asymptotic attractors for the 
dust era of evolution found above in eqn. (|l8|) providing the dust epoch has lasted long enough for the attractor to 
be reached. 

The behaviour of G, a and Ga in a universe like our own but which eventually becomes dominated by negative 
curvature is shown in Figure (0). 



VI. COSMOLOGICAL 'CONSTANT' ERA EVOLUTION 

In flat Brans-Dicke cosmologies a solution of the Friedmann equation with cosmic vacuum energy (p v = — p v ) is 

a = f BD+ % (30) 
<\> = <p o t 2 (31) 
Z2itp v 



(5 + 6ujbd)(3 + 2ujbd) 



(32) 



and p v is constant. This is not the general solution but it is the attractor for the general p v = —p v solution at late 
times J!?]]. It is a power- law inflation model Note that in Brans-Dicke theory, unlike in general relativity, a 

p v = —p v stress behaves differently in the Friedmann equation to an explicit constant A term [ |17| . It is the former 
that describes the stress contributed by a stationary scalar field with a constant potential. Every term in the BD 
Friedmann equations falls as t~ 2 for this solution. It is unusual in that it appears to predict that if the universe has 
just begun accelerating (as observations imply, then G should vary rapidly in the solar system. However, this 

argument assumes that the vacuum stress is dominant everywhere, right down to the solar system scale, which in 
reality it is not. 

If we substitute this solution for a(t) (but not 4>) in the ip and <j> evolution equations, (0) and (|5|) then we get, since 



that 



p v = —p v = const, 



2 , 3(2lobd + 1 ) 5 8tt -8tt ;a 

+ 2t * = 3TW 4P " ~ 2N ' ] * 3TW 4P " " ^ } 

■■ Z{2u BD + 1 ) ; 2 
V H 7T, "0 = e W CPm « 

It U) 



G 



So, at late times 



ip = Aa~ 3 = £)f-3(<"B D +i) 
ip = Er 3u)BD ~^ + F -> F 

3(2UJ BD + 1 ) • -87TCJ _g -327TW/9,, 

2t ^3T^ (4/, "- Qa } ^3T^' 



so 

<f> = A + Bt 2 + Cr 3wBD ~i -> O < 2 , (33) 

and, as expected, we get the same growing behaviour as in pure BD. When the universe becomes vacuum-energy 
dominated a tends to a constant value but aG cx G oc i~ 2 continues to fall. This behaviour is confirmed by numerical 
solutions shown in Figure (jl]). Using eqn. (|l8|), we see that if t v is the time when a vacuum-dominated era succeeds 
a sufficiently long dust-dominated era in a flat universe, then at t > t v in the vacuum-dominated era we expect 



a(t)Gf(*) S-f^ f . (34) 

Hence, today, we would have 

a(to)G(to) ~ CC " B J g ■ (35) 

7T0;(1 + Z v ) 1+2 "bd 

We see that, as in the situation where G is constant, the effect of a vacuum energy or quintessence field is to turn off 
variations in a when it takes over the expansion of the universe, H. 



VII. DISCUSSION 



We have formulated a simple gravity theory which extends general relativity, by the addition of two scalar fields, 
to include time variation of G and a. Previously, the study of the cosmological variation of physical 'constants' has 
confined attention to varying one constant only or to discussing the effects of altering the values of physical constants 
without a self-consistent theory for their dynamical variation p4l , |23| . The structure of unified gauge theories and 
particle physics theories with extra dimensions has given some indication as to the self consistency conditions required 
if traditional constants are allowed to vary, [^l],^2) . 

We have found that the expansion of the universe is affected by varying G to first order and the evolution of 
the expansion scale factor follows the behaviour found in Brans-Dicke cosmologies to leading order without being 
significantly affected by variations in a. The variations in a are affected by the variations in G through their influence 
on the expansion rate. This is significant in the dust-dominated era of cosmic expansion, which is known to exhibit 
a special mathematical behaviour in the absence of G variation, [Q. The effect of any G time variation simplifies 
the a variation and allows an exact solution to be found with a tx t n , where n = 2/(4 + 3ujbd) is determined by the 
Brans-Dicke parameter ujbd- 

In both the radiation and dust dominated eras, there are asymptotic solutions in which the product Ga remains 
constant and its value is determined uniquely by the coupling constants of the theory. However in universes like our 
own with the values of a and G near present values these asymptotic regimes are not reached throughout the life of 
the universe. Typically our present values for a are much larger than the values required by the asymptotic solution. 

In a curvature-dominated or quintessence-dominated era the variation in a ceases, just as in the situation with no 
G variation This is an important feature of all models with varying a in theories of the BSBM sort because it 
naturally reconciles evidence of variations in a at redshifts z ~ 1 — 3, with local (z = 0.1) constraints from the Oklo 
natural reactor if the universal expansion began to accelerate at z ~ 0.7, as current observations imply. 

Finally, we reiterate that the conclusions drawn above apply only to varying-a theories with negative £. The exact 
solutions given in eqs. jl5|),(|l6|),and ( p3| ) for the evolution of alt) during the radiation and dust eras no longer exist 
when C > and hence N < 0. During the curvature and cosmological constant-dominated eras the evolution of a{t) 
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is significantly changed by the variations of ip and the assumptions and ( pCf ) for the scale factor evolution are no 
longer valid. 

The study performed here provides a simple cosmological model in which the variation of two 'constants' can be 
studied exactly. A number of extensions are possible. The variations of weak and strong couplings can be included 
and the constraints imposed by any scheme grand unification can be imposed plj] . 
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